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NEOWISE Studies of Asteroids with Sloan Photometry: 

Preliminary Results 

A. Mainzer^ J. Masiero^ T. Grav^, J. Baucr^ ^ D. J. Tholen^, R. S. McMillan^ E. 
Wright^ T. Spahr^ R. M. Cutri^, R. Walker^ W. Mo^, J. Watkms^ E. Hand^ C. 

Maleszewski^ 

ABSTRACT 

We have combined the NEOWISE and Sloan Digital Sky Survey data to study 
the albedos of 24,353 asteroids with candidate taxonomic classifications derived 
using Sloan photometry. We find a wide range of moderate to high albedos for 
candidate S-typc asteroids that are analogous to the S-complex defined by pre- 
vious spectrophotometrically-based taxonomic systems. The candidate C-type 
asteroids, while generally very dark, have a tail of higher albedos that overlaps 
the S types. The albedo distribution for asteroids with a photometrically derived 
Q classification is extremely similar to those of the S types. Asteroids with similar 
colors to (4) Vesta have higher albedos than the S types, and most have orbital 
elements similar to known Vesta family members. Finally, we show that the 
relative refiectance at 3.4 and 4.6 /xm is higher for D-type asteroids and suggest 
that their red visible and near-infrared spectral slope extends out to these wave- 
lengths. Understanding the relationship between size, albedo, and taxonomic 
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classification is complicated by the fact that the objects with classifications were 
selected from the visible/near-infrared Sloan Moving Object Catalog, which is 
biased against fainter asteroids, including those with lower albedos. 



Introduction 



Asteroids and comets represent the leftovers from the formation of our Solar System. 
By studying their compositional variation, we can begin to better understand the conditions 
present at the earliest stages of planet formation as well as their subsequent evolution and 
processing. Asteroids are grouped into three main categories: C type or carbonaceous aster- 
oids, thought to be the most common type in the Main Belt; the S type or stony asteroids, 
a spectrally diverse group, and the X types, another diverse group of asteroids that have 
relatively featureless spectra but a wide range of albedos, probably representing a broad 



range of mineralogies and 
of the earlier taxonomic 


thermal histories. Gaffey et al. (1993 


) give an overview of some 


systems (e.g. Chapman, Morrison, & Zellner 


1975 Bowell et al. 


1978 Barucci et al. 1987 


Tedesco et al. 1989a Tedesco et al. 


1989b 


Howell, Merenyi & 


Lebofsky 1994). Tholen ( 


1984) created 14 taxonomic classes based on the combinations of 



colors available from the Eight-Color Asteroid Survey (ECAS; Zellner et al. 1985 Tholen 



& Barucci 1989). ECAS used filter passbands with wavelengths ranging from 0.34 to 1.04 
//m, since UV wavelengths were detectable by the photomultiplier available at the time. 
The Tholen and Tedesco systems also used albedo to differentiate asteroids; the X-group of 
asteroids have similar ECAS colors but markedly different albedos. 

The Tholen and Tedesco systems offer a powerful means of distinguishing unique asteroid 
groups. However, when CCDs replaced the photomultipliers at most observatories, they 
ushered in a factor of several improvement in quantum efficiency at red wavelengths compared 
to blue; while the ECAS photomultiphers extended out to 1.1 /xm, most CCDs' responses 
were diminished at this wavelength. These effects combined to make it difficult for observers 



to collect the full set of UV and near-infrared bands used by ECAS. Groups such as Xu et al. 



(1995), Lazzaro et al. (2004) and Bus & Binzel (2002a) obtained visible spectra of asteroids 



and used them for classification without albedos. Bus & Binzel (2002b) created a new 



taxonomic scheme relying solely on visible spectroscopy with 26 taxonomic types, and DeMeo 



et al. (2009) revised that scheme to 24 taxonomic types based on visible and near-infrared 



(VNIR) spectral signatures. Neither the [Bus fc Binzel] ( |2002b[ ) nor [DeMeo et al] ( |2009[ ) 
systems rely upon having albedo measurements for classification. Nevertheless, systems 
based solely on VNIR spectroscopy and photometry are widely used because observations at 
these wavelengths have generally been far more widely available than albedo measurements 
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for the ~500,000 asteroids known today. However, as discussed in Gaffey et al. (2002), the 



abihty to hnk taxonomic classification to asteroid mineralogy is complicated by the unknown 
surface properties such as particle size, the fact that some minerals have limited or no spectral 
features in the wavelengths used for classification, overlapping features, etc. It is therefore 
useful to try to understand how the various taxonomic types are linked to physical properties 
such as albedo and density. 

The fourth release of the Sloan Digital Sky Survey (SDSS) Moving Object Catalog 
(MOC; |Stoughton et al.|2"002||Abazajian et al.|2003| ) provided near-simultaneous observations 
of ~100,000 known asteroids in five bands {u, g, r, i and z) (Ivezic et al. 2001), and these 



have been used to study the distribution of colors throughout the Main Belt (c.f. Parker et al. 



2008 Nesvorny et al.|[2005 ). These studies have been conducted largely without reference to 
albedo, simply because the number of asteroids in the SDSS MOC vastly outnumbers those 



with well-measured albedos. To date, the Infrared Astronomical Satellite (IRAS; Tedesco et 
aL][2002| has been the largest source of radiometrically measured asteroid albedos, providing 
measurements for ~2000 asteroids. 



With the Wide-field Infrared Survey Explorer's NEOWISE project (Wright et al. 2010 



Mainzer et al. 2011a), thermal observations of > 157,000 asteroids are now in hand. In 



Mainzer et al. (201 Id), we compared the albedos derived from NEOWISE observations of 



~1900 asteroids with taxonomic types derived from VNIR spectroscopy. Here, we use the 
SDSS MOC photometry to obtain approximate taxonomic classifications for asteroids with 



NEOWISE observations. We initially focus on the system of Carvano et al. (2010 ), who define 



a classification algorithm based on SDSS colors that is compatible with previous taxonomic 
systems. 



2. Observations 



WISE surveyed the entire sky in four infrared wavelengths, 3.4, 4.6, 12 and 22 /xm 
(denoted Wl, W2, W3, and Wi respectively). Descriptions of the pre- launch mission design 



and testing can be found in Liu et al. (2008); Mainzer et al. (2005), and the post-launch 



description is given in Wright et al. (2010). A series of enhancements to the WISE data 



processing pipeline, known as NEOWISE, have enabled the detection of > 157,000 asteroids 
and comets throughout the Solar System, including the discovery of ~34,000 new minor 



planets (Mainzer et al. 2011a) 



24,353 objects, including nine NEOs and ~24,275 Main Belt asteroids (MBAs), were 
detected during the fully cryogenic portion of the NEOWISE survey and had matches with 
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SDSS MOC observations of sufficient quality to enable classification according to the method 
described in Carvano et al. (2010). The observations of these objects were extracted from 



the WISE archive using the First Pass version of the WISE data processing pipeline (Cutri 



et al.|[20TT ) following the methods and parameters given in Mainzer et al. (2011b hereafter 



M2), Mainzer et al. (2011c M3), and Mainzer et al. (2011d, M4). 



3. Preliminary Thermal Modeling 



We have created preliminary thermal models for each asteroid using the First-Pass Data 
Processing Pipeline described above. As described in references M2, M3 and M4, we employ 



the near-Earth asteroid thermal model (NEATM) of Harris (1998). The NEATM model uses 
the so-called beaming parameter t] to account for cases intermediate between zero thermal 
Lebofsky & Spencer (the Standard Thermal Model; 1989[ ) and high thermal inertia 



mertia 



Lebofsky et al. (the Fast Rotating Model; 1978); Veeder et al. (the Fast Rotating Model; 



1989); Lebofsky fc Spencer (the Fast Rotating Model; 1989). With NEATM, is a free 



parameter that can be fit when two or more infrared bands are available. Bands Wl and W2 
typically contain a mix of reflected sunlight and thermal emission. The flux from reflected 
sunlight was computed for each WISE band using the methods described in M2, M3 and 
M4; when sufficient reflected sunlight was present in bands Wl and W2, it was possible 
to compute the reflectivity at these wavelengths, piR, where we make the assumption that 
PiR = PsA = P4.6- The validity of this assumption and the meaning of piR is discussed in M4 
and will be the subject of future work. As described in M2 and M3, the minimum diameter 
error that can be achieved using WISE observations is ~ 10%, and the minimum relative 
albedo error is ~ 20% for objects with more than one WISE thermal band for which rj can be 
fitted. For objects with large amplitude lightcurves, poor H measurements, or poor signal 
to noise measurements in the WISE bands, the errors will be higher. 



3.1. High Albedo Objects 



We note that among the 24,353 asteroids considered here, there are ~55 that have 
Pv > 0.65. Of these, 48 have W3 peak-to-peak variations >0.3 mag, indicating that they are 
likely to be highly elongated. Almost all of the extremely high albedo objects have orbital 



elements consistent with membership in either the Vesta family or the Hungarias. Harris 



& Young ( 1988 ) and Harris et al. ( 1989 ) noted that E and V type asteroids can have slope 
values as high as G ~0.5. The assumption that we have used of G = 0.15 for an object like 
this would cause an error in the computed H for observations at 20° phase angle of ~0.3 
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magnitudes; this would drive the albedo derived using such an H value up by 0.3. These 
objects would greatly benefit from an improved determination of their H and G values. 



Discussion 



The Carvano scheme uses the SDSS colors as well as their measurement uncertainties 
to define nine spectral classes: Vp, Op, Sp, Ap, £„, Dp, Xp, Q^and C„. These are roughly 



analogous to the spectroscopically defined systems of Bus & Binzel (2002b) and DeMeo et 



al. (2009); the "p" indicates that the classification was derived photometrically. The nine 



Carvano classes were defined based on the ability of the SDSS colors to represent the system 



of Bus & Binzel (2002b). The exception to this is the Cp class, which is a conglomeration 



of the Bus K, L and Ld classes. For each SDSS observation, the probability that an object 
could be associated with a particular class was computed using the five SDSS magnitudes 
and their associated uncertainties; classifications and probabilities for ~63,000 asteroids were 



taken from Hasselmann et al. (2011). 



Figure [T] shows pv as a function of diameter for the individual Carvano classes. The 
bias of the VNIR SDSS survey against small, low albedo asteroids is evident in the dearth 
of objects in this regime, despite the fact that low albedo objects dominate most of the 



Main Belt (Masiero et al. 2011). As discussed in M4, this bias complicates study of the 



relationship between size and pv', this relationship is best studied using the full NEOWISE 
dataset as it is relatively unbiased with respect to pv- Figure |2] shows histograms of pv for 
the various taxonomic types, and Table 1 summarizes their statistical properties. The Sp 
types are fairly well-grouped but span a wide range of moderate to high py; while the Cp 
types tend to have significantly lower pv, roughly 15% have py>0.1, and the Dp types have 
36% with pv>0.1. In M4 we found that spectroscopically classified C types with diameters 
>30 km tended to have pv<0.1. The fact that so many Cp and Dp types with diameters 
<30 km have higher albedos suggests either that it is difficult to distinguish these objects 
with the SDSS colors in the Carvano scheme, that the bias against small objects with low 
Pv is skewing the result, or there is real scatter within the population of Cp and Dp types. 
The bias in the SDSS observations against small, low albedo objects means that measured 
Pv distributions will be skewed higher. Table 1 gives the statistical properties of the various 
classes both for all objects and for objects with classifications assigned with probabilities 
> 50%. As can be seen in the table, removing the objects with low probabilities does not 
significantly affect the results. 



Some asteroids have multiple observations in SDSS, and the different SDSS observations 
can each result in different classifications. Out of the 24,353 asteroids, 3924 received multiple 
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classifications. 



As described in Carvano et al. (2010), the most frequent variations occur 

As can be 



between classes with adjacent limits, such as Cp and Sp, and from Cp to Xp 
observed in Figure [2| the £Sp objects have an albedo distribution that resembles a blend 
of Sp and Cp types; similarly, the CXp types' albedos appear to be a cross between Cp and 
Xp types. Of more interest are the objects with classifications that span dissimilar classes, 
such as Cp to S'pOr Cp to Qp\ however, only a handful of objects in our sample fall into these 
classes. VNIR spectroscopy of these objects would illuminate what combination of slopes 
and/or absorption features gives rise to their diverse set of albedos. 



Cp and Dp Types. As discussed in M4, the B, C, D and T types of Bus & Binzel (2002b) 



and Tholen (1984) have similarly low pv, yet they are distinguished by markedly different 



Pir/Pv- We observe a similar phenomenon here with the Cp and Dp types. Figure |3] shows 
the distribution of pm/pv for those asteroids with sufficient reflected light in Wl or W2 to 
compute it. The Cp and Dp types have median py=0.064 ± 0.001 and 0.080 ± 0.002, re- 
spectively, yet the Cp types have p/R/py=1.147±0.032; Dp types have piji/pv='2.079±0.073. 
In M4, we hypothesized that the divergent Pm/pv was caused by fact that the Cp types' 
relatively fiat spectra from 0.4 to 1.0 /xm extends out to 3 — 4yum, while the Dp types slope 
steeply redward in VNIR wavelengths, continuing out to the Wl and W2 bands. Figure 
[5] illustrates the utility of pm/pv in distinguishing various "dark" asteroid types from one 
another. We have made the assumption that pwi = Pw2 at present and will revisit this 
assumption in the future. As discussed above, although the Cp and Dp asteroids generally 
have low albedos, they have a tail of higher albedos that are more similar to the bulk of the 
S-complex objects. Removing the objects with probabilities lower than 50% did not remove 
the tail of high-albedo Cp and Dp objects; this tail of high albedo objects could represent 
either real scatter within the populations' albedos or difficulty correctly distinguishing the 
C and D types on the basis of their SDSS photometry alone. 



X- Complex. 



The distribution of NEOWISE-derived albedos found for Xp asteroids 



closely resembles that observed in the general, undebiased Main Belt population (Masiero et 



al.|2011 ). We observed a similarly wide range of albedos in M4 for asteroids spectroscopically 
classified as X-complex. The wide range of albedos is not surprising, since the Tholen X 
type (from which the Bus and Bus-DeMeo X types are derived) consists of E, M and P 
asteroids that are separated by their albedos. The Carvano Xp type is based upon the Bus 
and Bus-DeMeo X types, neither of which use pv as part of their classification schemes. 



Sp, Cp, Ap, and Qp Types. Bus & Binzel (2002b) deemed the S-complex to have sufficient 



spectroscopic variation to warrant dividing it into six different subclasses. With the five SDSS 
bands, the Carvano system does not separate the S-complex. The Sp objects span a wide 
range of albedos with a similar though slightly higher mean value than the Cp types. Their 
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Pir/Pv values largely overlap, though the peak of the distribution shown in Figure[3]is slightly 
lower for the Cp types. The Cp type is an amalgamation of the Bus K, L and Ld classes, 
and these are distinguished by the degree of reddening shortward of i' band and flattening 
longward of ~1 ^m. The Ap types are thought to represent mantle material, and for the 
85 objects we observed with NEOWISE, we find that their albedos are very similar to the 
Sp types. Bus & Binzel (2002b) interpret Q-type asteroids as being the un-space- weathered 
parent bodies of the ordinary chondrite meteorites; see Chapman (2004) for a discussion of 
the effects of space weathering. In M4, we found NEOWISE observations of a handful of 
objects spectroscopically determined to be Q types, but their albedos were identical to the 
S complex; further, we found no significant differences between the DeMeo et al. (2009) S 
and Sw types, where the Sw type is thought to represent weathered S types. Here, we find 
424 Qp types with cryogenic NEOWISE observations sufficient to determine albedo, and our 
results reveal that the Qp types have slightly lower pv than the Sp types. All but two of the 
424 Qp types are Main Belt asteroids or Mars crossers, and they are distributed throughout 
the Main Belt as shown in Carvano et al. (2010); it would be worthwhile to obtain VNIR 
spectroscopy of these objects in order to determine whether or not the Q classification is 
accurate as the Main Belt is generally thought to be bereft of Q types ( Binzel et al.||2010 ). A 
nearly equal number of asteroids (453) are classified as being intermediate in type between 
Sp and Qpwith extremely similar py. 



Op Asteroids. Binzel et al. (1993) identified (3628) Boznemcova as the possible source 
of the L6 and LL6 ordinary chondrites based on the similarity of its reflectance spectrum to 
laboratory measurements of these common meteorites. Based on its spectral similarity to 
(1862) Apollo, they use Apollo's albedo and slope parameter G to compute a diameter of 7 
km. From NEOWISE observations of (3628) Boznemcova, we compute a diameter of 7.0±0.2 
km with an albedo of 0.330±0.043. Only two asteroids classified as O type spectroscopically 
according to the taxonomic system of Bus & Binzel (2002b) were observed by NEOWISE, 
(3628) Boznemcova and (169) Zelia, which is classified as an O type by Lazzaro et al. (2004) 
using the taxonomic system of Bus & Binzel (2002b). Using NEOWISE observations, Zelia 
has a diameter of 38.6±0.9 km and an albedo of 0.178±0.035; its ratio p//j/pv'=2.08±0.32. 



Carvano et al. (2010) identified 63 asteroids as 0-type candidates, and we have observed 16 



with NEOWISE, including Boznemcova. As shown in Table 1, these objects span a wide 
range of albedos, from 0.035 to 0.436, with a median value of 0.076; the largest is 12.1±0.2 
km in diameter, and the smallest is 1.75±0.23 km. None of these had measurements in 
Wl or W2, so Pir/pv could not be determined. It would be worth obtaining spectroscopic 
follow-up of these 0-type candidates (asteroids 11753, 20477, 21692, 29540, 32846, 34341, 
57736, 104596, 118595, 140949, 164099, 178653, 247450, 261454, 2005 VOlOl) to see if their 
spectra confirm their SDSS-based classifications. We note that the Op types with the highest 
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albedos also have the highest probabilities assigned by Hasselmann et al. (2011); all of the 
Op types with pv< 0-2 have probabilities lower than 25%. 



Vp Type. Carvano et al. (2010) found a clustering of Vp type asteroids near the Vesta 
family, with some scattered objects throughout the rest of the Main Belt. We observed 650 Vp 
objects with NEOWISE and find that they consistently have high py, tending to be brighter 
than the Sp types. Since the Vp types occupy a distinct region of albedo space compared to 
the background Main Belt asteroids (Masiero et al. 2011), it may be possible to use albedo 
in conjunction with orbital elements to find new candidate members of the Vesta family; this 
will be explored in a future work. The orbital elements of these objects are consistent with 
Vesta family members identified by Nesvorny (2010), with a median semi-major axis of 2.35 
AU and standard deviation (SD) of 0.13 AU; median eccentricity of 0.11 with SD=0.04; and 
median inclination of 6.5° with SD=2.4°. It would be desirable to obtain VNIR spectra of 
the Vp asteroids identified herein to further explore their true natures. 



Following the work of Ivezic et al. (2001 , 2002), we have plotted SDSS a* vs. i — z band 
magnitudes, with albedo shown as colored points (Figure [6]); the shading of Cp and Sp type 
points appears as expected from these prior works. As noted by Juric et al. (2002), asteroids 
with the bluest i — z colors appear to be associated with the Vesta family, and it can be seen 
in the figure that these objects tend to have high visible albedos. 



5. Conclusions 

By examining the intersection between the large thermal infrared NEOWISE dataset 
and the VNIR SDSS observations, we have further explored the relationship between tax- 
onomic classifications and pv for 24,353 asteroids and at 3-4 /xm for 1,819 asteroids. We 
find that the Carvano Sp and Vp types appear to be relatively robust in identifying aster- 
oids with high py. While the Cp types are generally significantly darker, we observe in 
the Cp and Dp populations high albedo tails that lead us to conclude that albedo cannot 
be conclusively determined purely from photometric VNIR taxonomy. We find a number 
of candidate 0-type asteroids among our sample; others have suggested that these are the 
progenitors of the LL6 ordinary chondrites, as well as Q-type candidates throughout the 
Main Belt. The objects with SDSS colors consistent with V-type asteroids typically have 
very high albedos, and their orbital elements are similar to known Vesta family members. 



suggesting that these new objects may indeed be family members as well. As in Mainzer 



et al. (2011d), we find that Dp types, which have a steeply red VNIR spectrum, also have 
systematically larger pm/py than Cp types, despite having similar py distributions. These 
higher pm/py values suggest that the Dp types' red spectrum continues through the WISE 
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Wl and W2 bands at 3-4 /xm. We caution that since these taxonomic types are determined 
from objects selected by a visible survey, the population is biased against low albedo objects, 
and hence the albedo distributions we have determined are similarly biased, particularly at 
the smallest size scales. This bias could be mitigated by obtaining a spectrally classified 
sample that better represents the true population of asteroids, both bright and dark. Such 
a sample could be created by choosing a set of objects from the NEOWISE survey, which 
has been found to be essentially unbiased with respect to pv (Mainzer et al. 2011e), and 
obtaining taxonomic classifications for these objects using a combination of literature values 
and new observations. While painstaking, such an effort would allow the distinction between 
observational bias and the physical links between size, albedo, and taxonomic classification 
to be clarified. 
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Fig. 1. — Visible albedo as a function of size for the nine individual classes defined by 
Carvano et al. (2010). The classes have been separated into two sub-panels for clarity. The 
bias of the SDSS survey against small, low albedo objects is evident. 
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Fig. 2. — NEOWISE-derived albedos of asteroids observed and classified by |Carvano et al. 



(2010). Only classes with more than ~100 asteroids are plotted, and these include some 
objects that have multiple SDSS observations that produced different classifications. The 
median py value is shown vertical red line. 
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Fig. 3. — NEOWISE-derived ratio Pir/pv ioi asteroids observed and classified according 
to the system of Carvano et al. (2010). Only asteroids for which Pm/pv could be fitted 
are included in this plot. While the Cp and Dp classes have similar pv values, pm/pv is 
distinctly different. The median pm/pv value is shown as a vertical red line. 
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Fig. 4. — NEOWISE-derived ratio Pm/pv vs. pv for asteroids observed and classified 
according to the Carvano taxonomic classification scheme. Only the ~1800 asteroids for 
which pm/pv could be fitted are included in this plot. 
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Fig. 5. — Orbital elements of the 650 asteroids classified as Vp] the typical ranges for the 
Vesta family are shown as dashed vertical lines (Nesvorny 2010). 
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Fig. 6. — SDSS a* vs. i — z colors for asteroids are shown with py and Pm/pv shown as 
different colored points. 



